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ABSTRACT 

We propose a new method for the characterization of stellar stratification in stellar systems. The method 
uses the mean-square radius (also called the Spitzer radius) of the system as a diagnostic tool. An estimate of 
the observable counterpart of this radius for stars of different magnitude ranges is used as the effective radius 
of each stellar species in a star cluster. We explore the dependence of these radii on magnitude as a possible 
indication of stellar stratification. This method is the first of its kind to use a dynamically stable radius, and 
though seemingly trivial it has never been applied before. We test the proposed method using model star 
clusters, which are constructed to be segregated on the basis of a Monte Carlo technique, and on Hubble Space 
Telescope observations of mass-segregated star clusters in order to explore the limitations of the method in 
relation to actual data. We conclude that the method performs efficiently in the detection of stellar stratification 
and its results do not depend on the data, provided that incompleteness has been accurately measured and the 
contamination by the field population has been thoroughly removed. Our diagnosis method is also independent 
of any model or theoretical prediction, in contrast to the 'classical' methods used so far for the detection of 
mass segregation. 

Subject headings: Galaxies: star clusters - Magellanic Clouds - stellar dynamics - Methods: statistical 



1. INTRODUCTION 

Mass segregation or stellar stratification is the phenomenon, 
according to which the massive stars in a star cluster are 
mostly concentrated toward its center. Due to dynamical evo- 
lution massive stars move inwards to the center of the cluster 
through two-body interactions. Due to the same process less 
massive stars may move outwards, so that on average the sys- 
tem 'evaporates' (e.g., Lightman & Shapiro 1978; Meylan & 
Heggie 1997). This is more commonly known as dynamical 
mass segregation. However, mass segregation is also being 
observed in dynamically young stellar systems, supposedly 
due to the formation of massive stars closer to the center of 
the newly-born system (e.g., Murray & Lin 1996; Bonnell & 
Davies 1998). This type of mass segregation is usually re- 
ferred to as primordial. 

Both types of stellar stratification affect the physical proper- 
ties of the cluster in a manner which allows us to develop diag- 
nostic tools for its detection and quantification. Such tools in- 
volve the identification of changes in the surface density pro- 
files and the mass functions of the clusters, or of the core radii 
of stars in different mass groups. Specifically, stellar stratifi- 
cation in a star cluster can be detected from the surface density 
profiles of stars in different brightness ranges. If / is the stel- 
lar surface density (for every brightness range) within a dis- 
tance R from the center of the cluster, then f(R) can be used 
as an indication of stratification at the magnitude limit, where 
it changes significantly (e.g., Scaria & Bappu 1981; Sagar et 
al. 1988;Subramaniametal. 1993; Kontizas et al. 1998). The 
density profiles can be approximated by log / oc 7 log R (El- 
son et al. 1987) and the slope, 7, can be used to quantify any 
significant change in f(R). Although this approximation was 
well designed, in practice peculiarities of the density profiles 
toward the center of clusters that are not completely spheri- 



cal tend to produce large uncertainties in 7, which hide any 
potential mass segregation (see, e.g., Gouliermis et al. 2004). 

Another well-applied method is the construction of the lu- 
minosity or mass function (LF, MF) of the cluster and the in- 
vestigation of its radial dependence (e.g., Pandey et al. 1992; 
King et al. 1995; Fischer et al. 1998; Kontizas et al. 1998, de 
Grijs et al. 2002a,b). If the cluster is indeed segregated, then 
the appearance of more massive stars toward its center will 
result in shallower LFs and MFs toward smaller radii (e.g., de 
Grijs et al. 2002a). However, this method has several con- 
straints imposed by the difficulties in the construction of an 
accurate MF based on the use of theoretical mass-luminosity 
(ML) relations. The use of different ML relations may result 
in significant differences in the MF slopes obtained, indepen- 
dently of whether any mass segregation is actually present (de 
Grijs et al. 2002b). In addition, the investigation of any radial 
dependence of the LF and MF slopes is not as straightforward 
as it seems, since its detection is rather sensitive to the radial 
distance range being considered, i.e., where the slope is to be 
estimated (Gouliermis et al. 2004). 

A rather interesting method is based on the use of the core 
radius of the cluster as estimated for stars in specific magni- 
tude (or mass) ranges (e.g., Brandl et al. 1996; de Grijs et 
al. 2002b). Core radii can be derived from the corresponding 
stellar number-density profiles by fitting to the latter the mod- 
els of Elson et al. (1987) and by approximating them to the 
canonical King models for Galactic globular clusters (King 
1966). However, this method seems to be model-dependent 
to a degree, which makes its application difficult. Another 
proposed method is to measure the mean stellar mass within 
a given radius (e.g., Bonnell & Davies 1998; Hillenbrand & 
Hartmann 1998), but this diagnostic requires accurate knowl- 
edge of the stellar masses, which is not always feasible (de 
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Grijs et al. 2002b; Gouliermis et al. 2004). 

For the methods discussed above, detecting stellar strati- 
fication requires observations at high resolution, so that de- 
tailed radial density distributions for the fainter magnitude 
ranges and complete MFs in several annuli around the center 
of the cluster can be constructed accurately. These methods 
can provide useful information about the segregated masses 
as well as on the radius up to which the cluster is shown to be 
segregated. However, it seems that all of these methods are 
very sensitive to the quality of the observed data, as well as 
to model assumptions. Moreover, all of the methods for the 
detection of stellar stratification used thus far are closely as- 
sociated to the fitting of predefined functions to the data and 
searching for parametric differences in these functions. One 
such method, for example, involves the study of the stellar 
projected density distribution of stars of different magnitudes; 
here, the differences in the exponents of the slopes fitted are 
seen as proof of mass segregation. A second method is based 
on the inherently imperfect derivation of radius-dependent 
MFs, while yet another makes use of the core radius of the 
system for different stellar groups, relying on apparent struc- 
tures, but it is extensively model-dependent. The problem 
that all these studies suffer from is that it is very difficult to 
compare their results among each other, thus leading to exten- 
sive discussions about the nature and even the reality of stellar 
stratification. 

In this paper we present a robust method for the detection 
of stellar stratification. The simple application of this method 
can yield direct information about the spatial extent (radial 
distribution) of stellar groups in different magnitude (mass) 
ranges. It is based on the notion of a dynamically stable ra- 
dius of a star cluster, the so-called 'Spitzer radius,' and leads 
to an observed 'effective radius.' The diagnostic method for 
stellar stratification is described in § 2. The diagnostic process 
is applied to and tested on a set of simulated spherical star 
clusters in § 3, which are constructed on the basis of Monte 
Carlo simulations. We thus validate the use of the diagnostic 
and the unavoidable constraints. We then apply the method 
to real data of two intermediate-age clusters in the Magel- 
lanic Clouds observed with the Advanced Camera for Surveys 
(ACS) on-board the Hubble Space Telescope (HST), and we 
report on our results in § 4. A comparison of the results of our 
method with those of the classical method of the MF depen- 
dence for the test clusters is also presented in § 4. Conclusions 
and a discussion on the use of this diagnostic tool are summa- 
rized in § 5. 

2. DIAGNOSTIC FOR SPATIAL STELLAR STRATIFICATION 

2.1. The Spitzer Radius 

The radius commonly known as the 'Spitzer Radius ' of a 
star cluster was introduced by L. Spitzer Jr. and his collabo- 
rators, who used it as a distance indicator for different stellar 
subclasses in a cluster (see, e.g., Spitzer 1958). Since gravity 
operates proportional to the inverse square of the distance be- 
tween a dynamical center and a subject mass, if the latter rep- 
resents a group of masses distributed around the center (e.g., 
stars in a cluster), then the square distance of this group can be 
replaced by the mean-square distance of its members. Specif- 
ically, the gravitational force in a spherical stellar system at a 
distance r from its center, is (e.g., Spitzer 1958, 1969) 



where V is the potential and M r the mass contained within a 
radius r. King (1965) showed that considering only one single 
type of stars with mass m and if S(x)dx is the number of stars 
in a strip of width dx, then assuming a Plummer (191 1) stellar 
distribution, the potential is 

F(r) = -— f S(x)dx. (2) 



Furthermore, Spitzer, investigating energy equipartition be- 
tween two different mass groups, notes that the term r 2 in 
Eq. (Q]) may be expressed as the mean value of the square dis- 
tance of all stars in a single mass group, if they were the only 
cluster members. Consequently, the square root of this value 
gives r in Eq. This radius, later called the 'Spitzer ra- 
dius,' corresponds to the distance up to which the stars of this 
specific mass group affect the gravitational field of the cluster 
as a whole. 

Therefore, the Spitzer radius of a star cluster is defined as 
the mean-square distance of the stars from its center, 
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where r-j is the radial distance of the i th star and N is the total 
number of stars. For a Maxwellian distribution of a group of 
masses in a parabolic potential well, which represents that of 
a globular cluster very well, the half-mass radius of the sys- 
tem is almost equal to (^90% of) the Spitzer radius (Spitzer 
1969). Consequently, considering the limitations in the mea- 
surement of the half-mass radius of a cluster versus the direct 
measurement of the Spitzer radius, the latter can be consid- 
ered quite important in relation to the dynamical status of the 
cluster, being used as its characteristic radius. 

As far as stellar stratification is concerned, in a star clus- 
ter displaying mass segregation one should be able to observe 
the more massive stars concentrated toward the center of the 
system by plotting the position coordinates of the stars ac- 
cording to the their magnitudes. The Spitzer radius, being a 
dynamically stable radius, can be used for the parameteriza- 
tion of the spatial distribution of stars in different brightness 
ranges. Based on this assumption, we propose here a method 
according to which stellar stratification in a star cluster can 
be quantified directly by estimating the corresponding Spitzer 
radius for every group of magnitudes, using Eq. (f3]). This sim- 
ple approach provides information on both where in the clus- 
ter every stellar-luminosity group is confined, and at which 
brightness stellar stratification occurs. 

2.2. Spitzer Radius and Stellar Stratification 

We test the hypothesis outlined above on artificial popu- 
lous spherical clusters that we constructed based on the use of 
the Monte Carlo technique. The proposed method is mainly 
designed for the identification of stellar stratification in clus- 
ters in the Magellanic Clouds (MCs), and therefore the simu- 
lated clusters are taken to have structural parameters similar to 
those in the MCs. Results based on observations of MCs clus- 
ters with HST are particularly considered due to their deep- 
ness and spatial resolution. Specifically, a sample of four star 
clusters (NGC 1818, NGC 2004, NGC 2100, and NGC 330) 
observed by Keller et al. (2000) with the Wide-Field Plane- 
tary Camera 2 (WFPC2) is used as an initial guideline in our 
simulations, by considering their structural parameters for the 
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artificial star clusters we construct here. As a consequence, 
our simulated clusters have tidal radii of rt ~ 2'.Q, while King 
models with C = log r t /r c ~ 1.0 are assumed to represent 
their density profiles. The corresponding core radii, r c , of 
the artificial clusters have values within the limits given by 
Mackey & Gilmore (2003a,b), who compiled two-color HST 
observations for a sample of 53 and 10 rich star clusters in 
the Large and Small Magellanic Cloud (LMC, SMC), respec- 
tively. 

The LFs of the simulated clusters are selected as varia- 
tions of a global LF, which was constructed based on the ob- 
served LFs of the clusters in the sample of Gouliermis et al. 
(2004). This global LF is found to be in very good agreement 
with the LFs of star clusters in the MCs as determined using 
HST imaging by other authors (e.g., de Grijs et al. 2002c: 
NGC 1805, NGC 1818, NGC 1831, NGC 1868, NGC 2209, 
and Hodge 14; Sirianni et al. 2002: NGC 330; Testa et 
al. 1999, Brocato et al. 2003: NGC 1866; Santiago et al. 
2001: NGC 1805, NGC 1818, NGC 1831, NGC 1868, and 
NGC 2209). An example of three artificial clusters with dif- 
ferent LFs is given in Fig.Q] The LF in each case is also given 
at the bottom panel of the figure. All three clusters have the 
same tidal radius and the same degree of stratification, but the 
differences in their LFs are fed through to their stellar num- 
bers and consequently to their appearance. 

Stellar stratification was considered for our artificial star 
clusters. Hence, several degrees of segregation (including no 
segregation) were adopted. This was done using a projected 
limiting radius within which each stellar group should be con- 
fined. Based on our hypothesis, stellar stratification in these 
clusters should be detected through the dependence of the es- 
timated Spitzer radius of stars in specific magnitude ranges 
versus the corresponding mean magnitude. Indeed, our tests 
show that the initially assumed stellar stratification (of any de- 
gree) is reconstructed by the plot of Spitzer radii per magni- 
tude range versus magnitude for all artificial star clusters con- 
sidered. In Fig. [2] we show a sample of three degrees of strat- 
ification for the same simulated cluster (similar to the middle 
panel in Fig.[T}. Spitzer radii in this figure are normalized for 
reasons of comparability. Three different symbols were used 
for the plot, as if the cluster was observed projected onto the 
xy, xz, or yz plane, respectively. 

2.3. Spitzer Radius of Observed Star Clusters 

Identical to the concept discussed in the previous sections, 
observations of the radial distribution of stars in a real star 
cluster give us a relation of the kind of Eq. ([3J, 
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However, this distribution is significantly affected by two im- 
portant observational constraints: (i) incompleteness of the 
stellar sample and (ii) contamination by field stars. To ob- 
tain an accurate measurement of the observed Spitzer radii per 
magnitude range in a cluster, both these constraints should be 
considered. 

2.3.1. Incompleteness of the observations 

The observed stellar samples in star clusters are incom- 
plete due to the observations. Although this can lead to a 
parameterization significantly differing from reality, the bias 



introduced to the data analysis is, in general, well under- 
stood. Consequently, solving the problem of incomplete- 
ness is a more or less straightforward procedure, which takes 
place through extensive artificial-star tests (e.g. de Grijs et al. 
2002a; Gouliermis et al. 2004). A set of artificial stars is gen- 
erated within each of the observed frames. Then, an identical 
reduction procedure is performed on the artificially enriched 
frames, in order to estimate the number of the 'recovered' ar- 
tificial stars. The completeness factor C is the ratio between 
the number of stars recovered to the number of stars origi- 
nally simulated. This factor depends on the brightness of the 
stars and their positions in the cluster. Therefore, an effec- 
tive completeness estimation for a stellar cluster requires that 
incompleteness is calculated for stars in different magnitude 
bins and for different distances from the center of the cluster. 

In order to correct the observed Spitzer radius for incom- 
pleteness, for every counted star its completeness factor is 
calculated according to its magnitude and radial distance. If 
the radial distance of the i th star is r, and its corresponding 
completeness factor is C, then we assign l/C, stars to this 
distance and the total number of counted stars is affected ac- 
cordingly. Thus, the completeness-corrected observed Spitzer 
radius should be 
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2.3.2. Contribution from the field populations. 

The contamination of the observations by the field popula- 
tions is probably the most difficult problem to be dealt with. 
Theoretically, if even one massive field star was counted as 
a cluster member and if it is far away from the center of 
the cluster, then the measured radius would be overestimated. 
In the same context, with comprehensive observations of the 
field we know the expected number of contaminating stars 
per magnitude range, but we still cannot know their position, 
which is the most important constraint as far as the radius es- 
timation is concerned. 

One may confront this problem by first defining the largest 
radial distance among the stars of the entire sample in a spe- 
cific magnitude bin and then normalizing the known surface 
density of field stars in this magnitude range to the surface 
confined by this radius. Hence, one could obtain a good es- 
timation of the number of field stars expected to contaminate 
this area in the specific magnitude range, and correct for field 
contribution by subtracting the appropriate fraction of the to- 
tal number of field stars. If, for example, there are Ny field 
stars expected in a total number, N, of stars in a specific mag- 



nitude range (within a maximum radial distance, 



<), then 



we can subtract their fraction F = Np/N < 1 from ev- 
ery counted star with distance r; < r max . Then, the field- 
subtracted observed Spitzer radius, ignoring for the moment 
any incompleteness corrections, derived from Eq. (0|i for stars 
in every magnitude range should be 
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FIG. 1 . — Sample of three different simulated star clusters using Monte Carlo simulations. The maps of the star clusters, which differ in their luminosity 
functions, are shown in the top panel. The corresponding LFs are shown in the bottom panel. All clusters were assumed to follow a King density profile and to 
have the same tidal radius of 0'.2. They are also similarly segregated. The LFs shown are in agreement with the results of various authors on star clusters in the 
MCs. 



However, this simplified approach assumes that the field 
stars are more or less homogeneously distributed within the 
boundaries of the cluster, which does not correspond to real- 
ity. Specifically, Eq. © weighs the contribution of each star 
to the sum of ri 1 by the probability, (1 — F), that it is a cluster 
member. Still, cluster members are by their very nature more 
centrally concentrated, and therefore stars at smaller radii are 
less likely to be field stars, while those at large radii have 
a higher probability that they are, in fact, field stars. As a 
consequence, in Eq. © true cluster members make a smaller 
contribution to the Spitzer radius than they should, while 
the contribution of the field stars is overestimated. We ver- 
ified this problem using additional Monte Carlo simulations 
for the construction of hypothetical background-field popula- 
tions. For these simulations, the use of a specific field LF was 
the only constraint considered. 

We constructed homogeneously distributed stellar popula- 
tions with LFs typical for the background fields of both Mag- 
ellanic Clouds, based on HST observations (e.g., Elson et al. 
1997; Holtzman et al. 1999; Dolphin et al. 2001; Smecker- 
Hane et al. 2002). Subsequently, we inserted our artificial 
star clusters (§ |2.21 i into the artificial background field. We 
constructed 32 different cases of four different clusters, each 
having four different degrees of stratification located in two 
different types of field (with different LFs). The application 
of our method to these artificial observations showed that it 
is indeed very difficult to disentangle the contribution of the 
field in the measured Spitzer radii on the basis of Eq. ||6). Con- 
sequently, and since it is not possible to quantify the effect 
of the field contribution to observed stellar samples in clus- 
ters as a function of distance from the cluster center, a more 



robust approach for the direct field decontamination of the 
observed stellar samples before measuring the Spitzer radii 
should be considered. Such a well-established method for the 
removal of any field contamination from the stellar samples 
of star clusters is based on the application of a Monte Carlo 
technique, which makes use of comprehensive observations 
of the local background field of the galaxy (e.g., Bonatto & 
Bica 2007). We apply such a sophisticated technique to actual 
HST/ ACS observations of two SMC clusters, for the estima- 
tion of their Spitzer radii as a function of stellar magnitude, in 
§3. 

2.4. The Effective Radius 

The method we propose for the diagnosis of stellar strati- 
fication in star clusters makes use of the observable counter- 
part of Spitzer radius, applied to the observed completeness- 
corrected stellar samples, and after the contamination of these 
samples by the field populations has been removed. This ra- 
dius, which we will refer to as the effective radius, r c s, is 
estimated for stars in specific narrow magnitude ranges; its 
dependence on magnitude is shown in the following sections. 
Eq. (0 can be reduced to the following expression and should, 
with all corrections properly applied, be equal to rg p . Thus, 
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where ?-j is the radial distance of the i th stellar member of the 
cluster (after field subtraction), in a specific magnitude range. 
Ci is the corresponding completeness factor at its distance for 
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FIG. 2. — Sample of three different degrees of segregation: (a) No stellar 
stratification, (b) smooth stratification, and (c) for more strongly segregated 
bright stars. The artificial cluster shown in the middle panel of Fig. |T]is used 
for this demonstration. The brightest (segregated) stars appeal' to be confined 
within smaller Spitzer radii than the faint (non-segregated) stars. Different 
symbols for each magnitude range represent rg p as it would be estimated if 
the cluster were observed projected onto the xy, xz, or yz plane, respectively. 



the same magnitude range and iV bs,cc,ff is me total number 
of observed stars corrected for incompleteness and reduced 
by the total number of field stars. 

This radius is a statistical quantity and its estimation is 
rather sensitive to the number of observed stars. Therefore, 
the effective radius can be easily over- or underestimated if 
this estimation is based on few stars only, which might be 
the case for the brightest magnitude ranges. In order to deal 
with this issue, one may apply the binning in variable magni- 
tude ranges and use wider magnitude bins toward the brighter 
limit, thus increasing the number of the sample stars. How- 
ever, with this solution one may not be able to observe the 
luminosity segregation of the stars in detail. It would be more 
appropriate to weigh each bin by the number of sample stars 
as an indicator of the statistical significance of the estimation 
of the corresponding effective radius. 

The uncertainty in the estimation of r e ff, as derived from 
Eq. ©, is mostly dependent on errors in the counting process, 
thus reflecting Poisson statistics. Another concern to be taken 
into account is the effect of projection on the estimated val- 
ues for the effective radius. For non-symmetrical loose stellar 
systems, it is found that projection does not significantly af- 
fect the true value of the Spitzer radius (e.g., Gouliermis et al. 
2000). Here, since we deal with spherically symmetric clus- 
ters, we consider the uncertainties introduced by projection 
for the simple case of a homogeneous density distribution. 

We assume a spherical distribution of N particles (N 3> 
1), uniformally distributed (with constant density, g) within a 
radius, R. If one observes this sphere, the distribution is not 
uniform anymore because of projection effects. Specifically, 
we investigate the size of the radius which includes half of the 
particles, projected along the line of sight, in relation to the 
radius of a sphere which - in reality - includes half of them 
(in 3D space). Because of symmetry we concentrate on one 
hemisphere, defined by the line of sight (direction of the Y 
axis in Fig.0. 

We define an elementary cylinder of thickness dx = 
R cos 9 dO, of base radius x = R sin (9, and of height 
y = R cos 9 (Fig. [3}. The angle 9 is measured in relation 




X 



YV 



FIG. 3. — Projection of a hemispherical distribution. The Y axis coincides 
with the line of sight. 



to the radius vertical to the surface of projection. The cylin- 
der contains dN = g dV particles, where V is the volume. 
Consequently, 



dN (9) = 2vr R 3 g sinfl cos 2 (9 d6. 



(8) 



Integrating over 0° < 9 < 90° we expect that the number 
of particles included is N = 2tt R 3 g/3. Half (N/2) are 
projected up to a limit cylinder base radius, x^, corresponding 
to f?h- This angle can be deduced from the relation 



N/2 = dN(6) = -2ttR 3 q / cos 2 6» d (cos 9) 



cos 3 # h 
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This angle gives = R sinf?h ~ 0.608i?. In reality, half 
of the particles are included in a volume contained within a 
sphere of radius R^ = 2~ 1 / 3 i? ~ 0.794i? (because of con- 
stant density). Thus, we have R^ = l.3xh. As a consequence, 
the proportion of the half-number radius (and consequently of 
the Spitzer radius, which is almost equal to the half-number 
radius for a relaxed stellar system) to its projected value for 
this distribution of particles is well-defined. 

In any case, in the present study this factor is irrelevant for 
further analysis, since the application of a constant conver- 
sion factor to the estimated effective radii is not necessary. 
Moreover, this application assumes that any cluster follows a 
homogeneous density distribution, which is of course not al- 
ways the case. Under these circumstances, we can treat the 
estimated effective radii used in the application of our method 
as the projected effective radii. 

3. THE STAR CLUSTER SAMPLE 

The proposed method for the detection of stellar stratifica- 
tion is designed primarily for the study of MC clusters, and 
therefore in the following sections we apply this method to 
such clusters to test its performance and to establish a consis- 
tent methodology. Moreover, in order to use the most com- 
plete available data we select clusters observed with the HST. 
The advantages HST introduced to crowded-field photome- 
try in the MCs have been documented extensively by sev- 
eral authors since its early observations and the quality of the 
data obtained with both the Wide-Field Planetary Camera 2 
(WFPC2) and the Advanced Camera for Surveys (ACS) has 
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FIG. 4. — Completeness functions for Kron 1 (solid line) and Lindsay 91 
(dashed line), respectively. The function has been integrated over position 
and color as a function of mp555w ■ 



proven to be more than adequate for detailed studies of young 
and old MC star clusters (e.g., Fischer et al. 1998; Johnson et 
al. 1999; Keller et al. 2000; Brocato et al. 2001 ; de Grijs et al. 
2002a,b,c; Stanghellini et al. 2003; Gouliermis et al. 2004; 
Mackey et al. 2006; Rochau et al. 2007; Xin et al. 2008). 

For the application of our new method for diagnosis of stel- 
lar stratification to real clusters we selected a set of data ob- 
tained with the ACS of the intermediate-age clusters Kron 1 
(Kron 1956) and Lindsay 91 (Lindsay 1958) in the SMC. Be- 
low, we discuss the reduction of the observations of these 
clusters (see also Mackey & Gilmore 2004), their photome- 
try, and the decontamination of the cluster populations from 
the contribution of the general background field of the galaxy 
using an advanced Monte Carlo technique. 

3.1. Data Reduction 

The SMC clusters Kron 1 and Lindsay 91 were observed 
with the Wide Field Channel (WFC) of ACS on-board the 
HST (program GO-9891, PI G. Gilmore). The ACS/WFC 
consists of two 2048 x 4096-pixel CCDs, separated by a gap 
of ~50 pixels. It covers a field of view of 202x202 arcsec 2 , 
with a scale of ~0.05 arcsec pixel -1 . The frames were taken 
in each of the F555 W and F8 14W filters. Exposure times were 
300 and 200 seconds, respectively. More details on the obser- 
vations and instruments can be found in Mackey & Gilmore 
(2004). 

The FITS files of the two clusters were retrieved from the 
Space Telescope Science Institute (STScI) Data Archive. The 
original data were reduced with the STScI pipeline, i.e., they 
had been bias and dark-current subtracted and divided by flat- 
field images. Photometry was performed with the ACS mod- 
ule of the package DOLPHOT 1 (Version 1.0). Photometry pro- 
cesses and corresponding parameters fully follow the proce- 

1 The ACS module of DOLPHOT is an adaptation of the photome- 
try pack age HSTphot (Dolp hin 2000). It can be downloaded from 

|http : / /pur cell ■ as ■ arizona . edu/dolphot /| 



dures and recommendation in the DOLPHOT manual. Photo- 
metric calibrations and transformations were done using the 
relations in Sirianni et al. (2005). 

We adopted three parameters in DOLPHOT to filter the 
photometric results, i.e., we selected only objects with 
—0.3 < sharpness < 0.3, crowding < 0.5 mag, and % 2 < 0.25 
in both frames. Meanwhile, we kept objects classified as good 
star (type 1) and star errors of types 1 to 7 by DOLPHOT, which 
are referred to as 'usable' in the DOLPHOT manual. To calcu- 
late the completeness of the photometry, DOLPHOT was run 
again in the artificial-star mode. For each cluster, we gener- 
ated ~ 10 6 fake stars with limits of 16.0 - 28.0 mag in bright- 
ness and —0.50 - 2.00 mag in color. The effect of unrealisti- 
cally large crowding that would lead to over-estimation of the 
incompleteness was taken into account by the DOLPHOT util- 
ity acsfakelist, which creates artificial star lists based on the 
original photometric catalog. The fake stars were binned in 
four dimensions, i.e, x and y positions, magnitude, and color. 
Fig. H] presents the completeness function of Kron 1 (solid 
line) and Lindsay 91 (dashed line), respectively. The func- 
tion is integrated over position and color and is a function of 
m F555W- More details on the data reduction are given in Xin 
et al. (2008). 

3.2. Field Decontamination 

The observed color-magnitude diagrams (CMDs) of the 
ACS fields centered on Kron 1 and Lindsay 91 are shown 
in Fig. [5] Both fields suffer from significant contamination 
by field stars, as can also be seen in the corresponding maps 
shown in Fig. [6] Without any information about cluster- 
membership probabilities for the observed stars, obtained, 
e.g., from radial velocities and/or proper motions, the quan- 
titative decontamination of the clusters from the field stars on 
a statistical basis becomes a fundamental method to obtain the 
CMD of the true stellar members of each cluster (see e.g. the 
pioneering work by Flower et al. 1980 on the LMC cluster 
NGC 1868). 

In this work, the algorithm described in Bonatto & Bica 
(2007) is used for field-star decontamination of the observed 
CMDs of both clusters. As a first step, we use the stellar 
number-density profiles of the clusters to identify the radial 
distance from the center of each cluster, Rii,„, where the stel- 
lar density becomes flat, indicating that the background-field 
density has been reached. We call this region the 'offset re- 
gion' of the cluster. All stars in the offset regions (with r > 
-Riim) are treated as field stars, while the rest (with r < i?n m ) 
are considered as the most probable cluster-member stars. 
Assuming a homogeneous field-star distribution, the number 
density of field stars is applied to the whole cluster region 
to remove the field contamination. We perform this calcula- 
tion in two dimensions in the CMD, i.e., in rr^ssw and in 
(niF555W — m F8i4w)> considering also the observational un- 
certainties in the photometry, <5p555w and <5fsi4W- 

In short, the decontamination process is done as follows: (i) 
we divide the CMDs of the cluster and offset-region stars, re- 
spectively, in 2D cells with the same axes along the mpsssw 
and (mp555w— niF8i4w) directions, (ii) we then calculate the 
expected number density of field stars in each cell in the CMD 
of the offset region, and (iii) we randomly subtract the ex- 
pected number of field stars from each cell from the CMD 
of the cluster region. In the following description we use the 
symbols x=m F 555w and £= (m F 555w-iriF8i4w) to simplify 
the notation. Consider a CMD cell with sides and coordinates 
X c ± Ax/2 and £ c ± A£/2, respectively, where (x c , £,c) are the 




cell's central coordinates. We assume a Gaussian probability 
distribution to calculate the probability of a star with CMD co- 
ordinates (X ± 6 X , £ ± 5 ( ), With ^=(^F555W +5 P814w) 1/2 > t0 
be found in the cell, and therefore the computations take into 
account both magnitude and color uncertainties, for example, 

P(x,x) = ^-e(- 1 / 2 )[(x-x)/^] 2 . (10 ) 
V 2tto x 

The expected field-star number density (pS? n ) in a cell is 
given by summing up the individual probabilities (P^ 11 ) of 
all offset-field stars (Nf s ) in the cell, divided by the offset area 
(A fc ), i.e., pf=Pf/A fs , where 

Jv fs „ . 

Pt l = J2l J p i(x,Xi^M)d X di. (ID 

i= 1 

-Pj(X) Xii £> &) I s me probability of the i field star, with 
CMD coordinates (x, and uncertainties (5 X i, 8^), to have 
the magnitude and color (x, 0- The integral is carried out in 
the two dimensions, \c- Ax/2 < X < Xc + Ax/2, and £ c : 
A£/2 < £ < £ c + A£/2, respectively. In the same manner we 
calculate the number density of the observed stars in the cell 
in the cluster-region CMD, p^s = ^"obs^cl, where A cl is the 
projected area of the cluster region (r < i?i; m )- Therefore, 
the expected number of field stars in the cell in the cluster- 
region CMD is n\f = (pg^/pgg) x n°g, where n^ 1 is 
the number of observed stars (at r < i?n m ) located in the 
cell. The number of probable cluster member stars in the cell 
will be ri™ 11 = — n^f 1 , and the total number of cluster 

member stars is the sum of all , i.e., N c \ = Ylcell n c\ U - 

To minimize any artificial effect intrinsic to the method, we 
apply the algorithm many times using different cell sizes (Ax, 
A£) and derive different decontamination results, from which 
we calculate the probability of a star to be identified as a 'true' 
cluster member. Finally the field-star-decontaminated CMD 



of the cluster contains the N c \ stars with the highest probabil- 
ity of being cluster members. Fig. [7] shows the correspond- 
ing CMDs of the cluster region (left panel), the offset region 
(central panel), and the final field-star-decontaminated CMD 
(right panel) for both clusters. 

4. APPLICATION OF THE DIAGNOSTIC METHOD 

The proposed method for the diagnosis of stellar stratifica- 
tion is based on the assumption that the effective radius of 
stars in a specific magnitude (or mass) range will be a unique 
function of this range if the system is segregated. In this sec- 
tion we apply our new method for the detection and quantifi- 
cation of stellar stratification to the intermediate-age clusters 
Kron 1 and Lindsay 91. Before the application of our diag- 
nostic method and to check if the clusters are segregated, we 
investigate the stellar stratification using one of the most accu- 
rate 'classical' methods, the radial dependence of the cluster 
LFs and MFs. 

In the previous section we showed that both clusters are 
evolved, by virtue of the presence of clear red-giant branches 
(RGBs) and red clumps in their CMDs. However, in the cor- 
responding 'clean' CMDs, obtained after decontamination of 
the cluster stellar samples by the field populations, shown in 
FigE] (right panels), it can be seen that there are stars at the 
upper main sequence (MS) and just above the MS turn-off, 
which are not sufficiently removed from the observed CMDs. 
These stars are most probably members of the younger SMC 
field, to which both clusters belong, and they will not be con- 
sidered in the following analysis. As a consequence, for the 
investigation presented here, we use only the lower-MS popu- 
lation below the turn-off of each cluster and that of the RGB. 

4.1. Stellar Stratification in the selected clusters 

In this section we apply the most effective classical diag- 
nostic method developed in the past for the investigation of 
stellar stratification to our ACS photometry of Kron 1 and 
Lindsay 91, to check if the clusters are indeed segregated. 
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FIG. 6. — The observed ACS/WFC fields of view of Kron 1 (left) and Lindsay 91 (right). The overlayed circles have radii, Rn m , corresponding to the boundary 
between the clusters and the gene ral f ield. They were defined based on the clusters' surface density profiles, as the distance from the cluster center where the 
density profiles become flat (see § 13.21 . 



This well-applied method requires the construction of the 
cluster LF and/or the MF and the investigation of its radial de- 
pendence. If the cluster is indeed segregated then the appear- 
ance of more massive stars toward its center will result in a 
shallower MF at smaller radii. De Grijs et al. (2002b) present 
the difficulties in the construction of an accurate MF based on 
the use of different ML relations and stress the differences in 
the MF slopes obtained. Moreover, significant numerical bi- 
ases in the determination of the slope of the MF using linear 
regression are found if the construction of the MF was made 
from uniformly binned data (Mafz Apellaniz & Ubeda 2005), 
implying the presence of systematic errors in the slopes of 
MFs calculated in this way. 

Gouliermis et al. (2004) argue that the investigation of any 
radial dependence of the MF slope is not as straightforward as 
it seems, since any such dependence is rather sensitive to the 
selection of the radial distances being considered. Although 
stellar stratification can be identified from the cluster LF, its 
quantification is quite difficult because the approximation of 
the LF with a single power law, which is the usual approach, 
can only give very rough results. In any case, we apply this 
diagnostic to our data, in order to identify stellar stratification 
in our clusters and to be able to compare the results of a well- 
established classical method with those of our method. 

4.1.1. Radial Dependence ofLFs and MFs of the Clusters 

In this section we examine the dependence of the shape and 
slope of the stellar LF on position within the clusters. We ba- 
sically use the cluster LFs as smoothing functions of the full 
2D (mF555w, "1F555W — 7™F8i4w) CMDs shown in the right- 
hand panels of Fig. [7] In Fig. [8] we show all annular cluster 
LFs out to r = 100.0", corrected for the effects of incom- 
pleteness and background contamination at each respective 
annular radius. We carefully considered the optimal radial 
ranges to be used for our LF analysis, which was in essence 
driven by the need to have statistically significant and simi- 
lar (cf. Mafz Apellaniz & Ubeda 2005) numbers of stars in 



each of the subsamples used for our comparison of the differ- 
ent annular LFs. For reasons of clarity, we did not normalize 
the LFs to the same sampling area since this would result in 
smaller separations among the various curves, hence hamper- 
ing our assessment of any radial dependence of the LF shape 
(at this point we are not interested in the absolute numbers of 
stars per unit area). 

We note that for a proper comparison of the annular LFs 
in Fig. [8] one should only consider the magnitude range in 
between the vertical dotted lines. For brighter stars, stochas- 
ticity in the stellar LFs starts to dominate, while the LFs are 
also significantly affected by the presence of red-clump stars, 
which are clearly visible in both CMDs. For fainter stars 
our subsamples are significantly statistically incomplete, at 
<50%. On the basis of a casual comparison of the annular 
LFs in Fig. [8j we conclude that there is little evidence for 
luminosity segregation in Kron 1, although there is a hint of 
a flattening of the LFs towards smaller radii, particularly for 
WF555W ~ 22.5 mag. Lindsay 91, on the other hand, ex- 
hibits clear signs of luminosity segregation for most of the 
magnitude range of interest. 

The conversion of an observational LF in a given passband 
i, <p(Mi), to its associated MF, £(m), is not as straightforward 
as often assumed (see, e.g., de Grijs et al. 2002b). The dif- 
ferential present-day stellar LF, dN/d<p(Mi), i.e. the number 
of stars in the absolute-magnitude interval [Mi, Mi + dMi], 
and the differential present-day MF, diV/d£(m), i.e. the mass 
in the corresponding mass interval [in, m + dm], are related 
through dA = —<p(Mi)dMi = £(m)dm, and therefore 

<t>{Mi) = -^H^- ( 12 ) 

Thus, to convert an observational LF into a reliable MF, one 
needs to have an accurate knowledge of the appropriate ML 
- or mass-absolute-magnitude - relation, dm/dMj (for a de- 
tailed discussion, see de Grijs et al. 2002b). In fact, it is the 
slope of the ML relation at a given absolute magnitude that 
determines the corresponding mass, which is therefore quite 




FIG. 7. — CMDs of the cluster regions (r < i?n m ; left panel), the offset regions (r > i?n m ; central panel), and the final CMDs of the clusters after field-star 
decontamination (right panel) for both Kron 1 (top) and Lindsay 91 (bottom). 



model dependent. This has been addressed in detail by, e.g., 
D'Antona & Mazzitelli (1983), Kroupa et al. (1990, 1993), 
Elson et al. (1995) and Kroupa & Tout (1997). Given the non- 
linear shape of the ML relation (de Grijs et al. 2002b) and the 
small slope at the low-mass end, any attempt to model the ML 
relation by either a polynomial fit or a power-law dependence 
will yield intrinsically unreliable MFs (cf. Elson et al. 1995, 
Chabrier & Mera 1997), in particular in the low-mass regime. 
This model dependence is clearly illustrated by, e.g., Ferraro 
et al. (1997) and de Grijs et al. (2002b), who compared the 
MFs for their sample clusters derived from a variety of differ- 
ent ML relations at that time available in the literature. Nev- 
ertheless, for the sake of our discussion on mass segregation 
in our two sample clusters, we only need to consistently apply 
the same ML relation to the clusters' annular LFs; any differ- 
ences in the resulting MF shapes will then be due to intrinsic 
differences in the stellar mass distributions as a function of 
radius in the clusters. 

In addition, the true structure of the CMD above the MS 
turn-off is very complex. Stellar populations of different 
masses overlap in color-magnitude space, so that unambigu- 
ous mass determination from isochrone fits, for the few dozen 
stars populating these areas in each cluster, is highly model- 
dependent. On the other hand, while different evolutionary 
models differ in terms of absolute calibration of the mass- 
luminosity relation, they show a quite good agreement in 
terms of relative mass distribution in well established evolu- 



tionary phases. In any case, in a differential analysis such 
as presented here, the uncertainties involved in the mass de- 
terminations of the evolved stars are too large and system- 
atic, so that we cannot include these stars in our MF anal- 
ysis. Therefore, we used all stars below the MS down to 
the 50% completeness limit of the data. After applying 
the appropriate distance modulus to obtain absolute magni- 
tudes, for our luminosity-to-mass conversions we used the 
updated Padova isochrones (Marigo et al. 2008; available 
from |http : //stev. oapd. inaf . it/cmdj l for the appropri- 
ate HST photometric system, the metallicity and age of the 
cluster, and a Kroupa (1998) IMF corrected for the presence 
of a binary population, as suggested for use with the web in- 
terface. For the metallicity and age of both Kron 1 and Lind- 
say 91 we used Z — 0.001 (where Z = 0.020) and an age 
of 6.3 Gyr; their distance moduli are m — M = 19.34 mag 
and to — M = 19.10 mag, respectively. These physical pa- 
rameters are derived from fitting the isochrone model that best 
matches the observed features of the CMD of each cluster. 

The MS turn-off magnitude in both Kron 1 (TOF555W — 
22.3 mag) and Lindsay 91 (tofsssw — 22.1 mag) corre- 
sponds to a mass of log^^/M©) ~ —0.009, while the 50% 
completeness limits correspond to log(TO*/M ) = —0.2 and 
—0.17 for these clusters, respectively. In Fig. [9] we show 
the derived MF slopes as a function of cluster radius for our 
adopted ML conversion in the full mass range; the canonical 
Salpeter slope would be a — —2.35 in this representation. 



10 



J_j 1-1 

a 

a 



as 



Gouliermis, de Grijs, & Xin 

2.6 



Kron 1 



s - 



stochastic affacbs 
dcUilUL* (HC) 



jit 

! /*'// 

■ :/ */ H 

n i '/ t< 

it tl/f /.' 

, t ft 

It Si ,'- I " 



1 3 

i 

i 



> \ VI • 



», 

I 

I 



black: r 

nil r 

blue: r 

green; r 

eyan: r 



[o.so]" 
■■i.-i 
[40.60]" 

[ea.BDi" 

[80. 1 OS]" 



£Mol normalised by sts») 



£8 

m F5(55W ( ma g) 



S 



J, 

a; 

130 

o 



0.6 



i 1 1 1 1 1 1 1 1 1 [ 1 1 — r-> r 

Lindsay 91 



stochastic effects : 
dominate f,RC) 



1 1 



1 jv 



Hi 



/■ 
( 



alack: r = [0.30]" 
red: r - [B0.+0)" 
Bills: r - [40.BO]" 
green: r = [S0.60]" 
cyan; r = [80.1O0]" 
(Not normaliaad by area) 



20 



24 



(mag) 



as 



FIG. 8. — Completeness-corrected LFs of Kron 1 (left) and Lindsay 91 (right) within different annuli around the center of the clusters. These LFs clearly show 
how difficult it is to derive a single slope for the entire LF to compare the results for different clusters. Specific parts of the LF can, however, be approximated by 
a single power law. The LF of Kron 1 does not show any prominent dependence on the radial distance. Lindsay 91, however, shows a clear radial dependence of 
its LF in the magnitude range 22.0 ^5 mp555w/rnag 25.0. One should only consider the LFs in the magnitude range in between the vertical dotted lines, as 
the brightest stars are affected by small-number statistics and the presence of a red clump (RC), while the lowest-luminosity stars suffer from significant statistical 
incompleteness. 



The adopted radial ranges for our annular MFs are indicated 
by the horizontal 'error' bars. The vertical error bars represent 
the formal uncertainty in the fits. Although the uncertainties 
are large, both clusters seem to be affected by mass segrega- 
tion, with Lindsay 91 being most obviously mass segregated, 
out to at least r ~ 80". 

4.2. Effective Radii of the Clusters 

We calculated the effective radii of the stars in different 
magnitude ranges for both clusters. Magnitude bins of vari- 
ous sizes were tested and we found that a reasonable bin size, 
which provides good statistics and allows for the detailed ob- 
servation of the dependence of r e g on magnitude, is 0.5 mag. 
We calculated the effective radii of the stars as a function of 
magnitude in the F8 14W filter. The relations derived between 
the computed effective radius and the corresponding magni- 
tude bin for both Kron 1 and Lindsay 91 are shown in Fig. 
[TOl It is shown that, in general, the effective radius behaves 
as a function of magnitude for both clusters, thus providing 
clear indications of stellar stratification. While this behavior 
is more definite for the fainter stars, where it is seen that for 
both clusters the effective radii of these stars are larger, the 
uncertainties in the r e g calculation for the brightest stars are 
rather large due to small-number statistics, and therefore no 
definite trend for the bright stars can be derived. 

From the r e g vs. magnitude plots in Fig.[l0]one may con- 
clude that the clusters exhibit stellar stratification but to a dif- 
ferent degree. Moreover, from the plots in Fig. QJJone can 
derive the degree of stratification of the clusters in terms of 
the brightness range of the segregated stars, as well as the ef- 
fective radius within which they are confined. Considering 
that a steeper slope of the relation r c fj (tofsi4w) represents a 
higher degree of stratification, we conclude that Lindsay 91 is 



more strongly segregated than Kron 1 . 

If we want to parameterize the differences between the clus- 
ters and obtain more quantitative results, we should use the 
two primary output parameters of our method, i.e., the mag- 
nitude and the effective radius of segregation. Both represent 
the limit (in magnitude range and radius, respectively) beyond 
which the slope of r e fi(mF8i4w) changes significantly. How- 
ever, as can be seen in Fig. \W[ r c g is not a mono tonic func- 
tion of brightness. Specifically, for Kron 1, while for stars 
with TOp8i4W r$ 19 mag there is a trend of r e g to become 
larger for fainter magnitudes, this trend cannot be confirmed 
statistically due to large uncertainties. For fainter magnitudes, 
down to ~ 20.5 mag, the relation r e g vs. TOF814W shows 
fluctuations and no specific trend. For even fainter stars, with 
TOF814W ~ 21 mag, the slope of this relation shows a definite 
steepening as a clear indication that these stars are indeed seg- 
regated. However, this slope is still quite shallow, providing 
evidence of a low degree of segregation. A comparison be- 
tween the r e g vs. tofsmw relation of Fig.[l0]for Kron 1 and 
the LFs of the cluster constructed for different distances from 
its center (Fig. [8), shows that r e s seems to be a function of 
brightness for the entire observed magnitude range, but statis- 
tics do not allow the verification of a definite relation. Specif- 
ically, Fig. [8] shows that the LF slope seems to be distance- 
dependent for the full extent of the cluster and for stars within 
the entire observed brightness range. Indeed, a functional re- 
lation between r e g and brightness seems to exist in Fig. [10] 
for the entire cluster, with r e g ~ 0.7' for mF8i4W ~ 17 mag, 
but this trend is not supported statistically. The small numbers 
of the brightest stars allow a solid statistical interpretation of 
this plot only for stars of mpsiAW ~ 21 mag. For these stars 
the LF slope becomes steeper outwards. The advantage of 
our method lies in the fact that we are able to define in a di- 
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FIG. 9. — Radial dependence of the slopes of the differential MFs of Kron 1 
(top) and Lindsay 91 (bottom). This dependence provides a clear indication 
of mass segregation of the stars on the lower MS (below the turn-off) for both 
clusters, and particularly for Lindsay 91. However, the uncertainties in the 
counting process, the model dependence for construction of the MF, and the 
selection of specific radial distances whithin which the MF is constructed for 
the analysis of stellar stratification introduce additional uncertainties, which 
work their way through to the results of this analysis. 



rect manner the distance of segregation for these stars. This 
distance, derived from Fig.[l0l is r e g ~ 0.8'. 

For Lindsay 91a more obvious trend of r e g as a function 
of brightness can be seen for almost the entire observed mag- 
nitude range, mpsuw ~ 18.5 mag. Also in this cluster, no 
clear correlation between r e g and mpsi4w is observed for the 
brighter stars due to their small numbers. The relation be- 
tween r c ff and mF8i4W can be clarified easily, however, for 
stars with mpgi^vv ^ 21 mag. It is interesting to note that 
from the comparison of the relation of r c s vs. mp8i4w of 
Fig.[l0]for Lindsay 91 with the corresponding LFs for differ- 
ent distances from the cluster center, shown in Fig. [8] we see 
that both plots agree that significant segregation is observed 
for stars with tofsi4W ~ 21 mag. However, specific radial 
distances were selected for the LFs of Fig. [8] and therefore 
we cannot define an accurate distance of segregation for these 



stars, except that segregation occurs at distances r ?J 40". 
With our new method we can accurately define the distance 
of segregation. From Fig.[l0]we derive that this distance cor- 
responds to the effective radius of r c g ~ 0.72' = 43", which 
agrees very well with the results from the LFs method, but 
only because of the arbitrary selection of this distance for this 
method. 

It should be noted that the relation of r c g vs. magnitude 
derived for the MS stars (below the turn-off) is the most reli- 
able indicator of the degree of segregation of the studied clus- 
ters. In the case of this relation for the brighter evolved stars 
(above the turn-off) the situation is more complicated because 
the mass of these stars changes very slowly with magnitude. 
Specifically, according to the evolutionary models adopted 
here (Marigo et al. 2008) and the isochrone of 6.3 Gyr, the 
difference in mass of stars evolving from the base of the RGB 
to the AGB is much less than 0.01 M©, and thus the corre- 
sponding values of r ff refer to extremely narrow ranges of 
mass per magnitude bin. As a consequence different magni- 
tude bins correspond to roughly the same mass and the rela- 
tion of r c ff vs. magnitude for these stars provides indications 
of stratification of of stars of different magnitude but nearly 
identical mass. Naturally, this effect has consequences on the 
detection of mass segregation. In the case of the relation of 
r e ff vs. mass, the fact that the evolved stars have small mass 
differences can improve the number statistics, since a larger 
number of stars will correspond to the same mass-bins. In- 
deed, considering that according to the isochrone model all 
stars brighter than the turn-off have masses around 1 M©, if 
we group together all these stars, then we should combine the 
first 10 magnitude bins shown in the plots of Fig.[l0]into one. 
The corresponding values of r c g are 0.80' ± 0.06' for Kron 1 
and 0.88' ± 0.07' for Lindsay 91; both values are quite con- 
sistent with the derived trends. 

In general, based on the above, with the method for detect- 
ing and characterizing the effect of stellar stratification pro- 
posed here we can distinguish between different degrees of 
segregation among star clusters and compare our results from 
one cluster to another. The values of r c ff derived by using only 
MS stars (below the turn-off) for the studied clusters should 
be considered as the best tracers of mass segregation. 

5. CONCLUSIONS 

We present a robust new method for the diagnosis of stellar 
stratification in star clusters. This method uses the effective 
radius, which represents the observed Spitzer radius as esti- 
mated for various stellar groups of different magnitude (or 
mass) ranges in the same system. This radius is an impor- 
tant dynamical parameter, which approximately describes the 
area of gravitational influence of a given stellar group in the 
cluster. The proposed method detects stellar stratification in 
star clusters: if a cluster is segregated, then the effective radii 
for different stellar groups should be different from one group 
to another. A relation between the effective radius of every 
group and the corresponding selected magnitude (or mass) 
range can be established to consistently characterize any strat- 
ification. As a result, comparison among the results obtained 
for different clusters can be made. 

Stellar stratification in star clusters is usually quantified 
based on either the surface density profiles of stars in differ- 
ent magnitude ranges or the cluster LF (or MF) at different ra- 
dial distances from its center. Using the first method one may 
observe the magnitude (mass) limit for the segregated stars 
and using the second the radial distance where they are ob- 
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FIG. 10. — Estimated effective radii of stars within different magnitude ranges versus the corresponding mean magnitude for both clusters. 



served to be segregated. Both of these results can be derived 
using the single application of the method proposed here. Fur- 
thermore, to apply the 'classical' methods for the detection of 
mass segregation, certain complicated fitting procedures must 
be applied, leading to model-dependent results. On the other 
hand, our method is much more straightforward in its applica- 
tion, and more precise in its results. This is mostly so because 
it is not necessary to (i) introduce any significant assumptions 
to estimate the cluster parameters, and (ii) apply any fit to 
parameter correlations and subsequently test the relation be- 
tween the model output and the observables. Both of these 
steps must be applied in the classical approaches to character- 
ize stellar stratification. 

We present the application of this new diagnostic tool 
to two mass-segregated clusters in the SMC, observed with 
HSTI ACS. The application of the new method to space-based 
observations helped us to establish a scheme for the compar- 
ison between the results for different clusters and to define 
the important parameters required to characterize stratifica- 
tion. These parameters are the magnitude and radius of seg- 
regation, which specify the limits beyond which segregation 
becomes stronger and more apparent. In addition, one may 
quantify the differences between clusters through the slope of 
the relation between the effective radius and the correspond- 
ing magnitude range, which can be used as an indicator of the 
degree of segregation of each cluster. 

In each of the clusters observed with the HST, considered 
here, it was found that segregation is apparent for stars over 
the full observed brightness range, but small-number statis- 
tics do not allow the detection of a definite trend between 
r e g and magnitude for several of the brightest stellar groups. 
A comparable radius of segregation is observed for stars of 



wiF8i4W <; 21 mag in both clusters, but the degree of segre- 
gation seems to differ from one cluster to the other, with Lind- 
say 91 exhibiting a steeper function of r e ff(mF8i4w) than 
Kron 1. 

The continuous increase of r e g as a function of brightness 
away from the radius of segregation up to the observed lim- 
its of the clusters shows that the fainter stars are distributed 
throughout a larger volume, an observation reminiscent of the 
'evaporation' of star clusters. Still, these results are based on 
only two clusters and can therefore only be used rather quali- 
tatively. 

We conclude that the method we present here is simpler 
and more accurate than previously developed methods. This 
is based, first and foremost, on the more straightforward way 
in which the effective radius is estimated and, secondly, on 
avoiding any model-dependent fits. Instead, a direct correla- 
tion between the cluster's observables is sufficient to exhibit 
any stellar stratification and to parameterize its significance. 
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